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Abstract 

Using the theory of combinatorial species, we compute the cycle index for bipartite graphs, which we 
use to count unlabeled bipartite graphs and bipartite blocks. 



1. Introduction 



o 
u 

1 -p^ . A hicolored graph is a graph each vertex of which has been assigned one of two colors such that each edge 

C^ ' connects vertices of different colors. A bipartite graph is a graph G which admits such a coloring. Given 

j white and k black vertices, there are 2^^ ways to join vertices of different colors. Thus the number of 
(labeled) bicolored graphs on n vertices is 

>: &„- V i%'^. (1) 



i+j— n 



Bipartite graphs are not so easy to count directly. Every connected bicolored graph has exactly two colorings 
in white and black, so we can count bipartite graphs by relating them to connected bipartite graphs. To 
■"sj" , do this, we use the exponential formula [T^ section 5.1], which implies that if f{x) — X^riLi /n^;"/'*' i^ 

^^ ' the exponential generating function for a class of connected graphs then e^^''^' is the exponential generating 

function for graphs all of whose connected components belong to this class. It follows that with B{x) = 
T^'^=Q hnx" /nl, where &„ is given by equation yj), the exponential generating function for connected bicolored 
graphs is logi?(x), the exponential generating function for connected bipartite graphs is i logi?(a;), and the 
^^ I exponential generating function for bipartite graphs is e'°^^^^-'/^ — y^B{x). 

^ ■ Just as arbitrary graphs may be decomposed into their connected components, arbitrary connected 

graphs may be decomposed into "blocks" — maximal 2-connected (or "nonseparable") subgraphs. Techniques 
developed by Robinson [13] were applied by Harary and Robinson [9] to show that the exponential generating 
function N{x) for labeled 2-connected bipartite graphs is related to the exponential generating function 
P(x) = ilog_B(a;) for connected bipartite graphs by the equation logP'{x) = N'{xP'{x)). This equation 
suffices to compute the number of labeled bipartite blocks on n vertices and their asymptotics. 

To count unlabeled bipartite graphs we can take a similar approach. It is not too difficult to find the 
generating function for bicolored graphs from first principles, in a way that is very similar to counting 
unlabeled graphs (see, e.g., [7]). There is an analogue of the exponential formula for unlabeled graphs (see, 
e.g., [HI equation (3.1.1)] and [H p. 46, equation (20b) and p. 55, equation (60 ii)] so we can easily relate the 
generating function for bicolored graphs to that for connected bicolored graphs and the generating function 
for connected bipartite graphs to that for all bipartite graphs. The difficult step is relating connected 
bicolored graphs to connected bipartite graphs: some unlabeled connected bipartite graphs can be bicolored 
in two different ways, and some in only one way, as shown in Fig. [l|. 
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(a) A connected bipartite graph with two distinct (b) A connected bipartite 

bicolorings, one of which is shown graph with just one distinct 

bicoloring 

Figure 1: Some connected unlabeled bipartite graphs have two distinct bicolorings, but some have only one 

To deal with this problem, we consider the two-element group 62 acting on bicolored graphs by inter- 
changing the colors. We want to count orbits of this group acting on connected bicolored graphs. To do 
this using Burnside's lemma, we would need to know the number of connected bicolored graphs fixed by 
each of the elements of ©2- This is not so easy to compute directly; however, it is not hard to compute the 
corresponding information for the action of S2 on all bicolored graphs, and by using an extension of the 
exponential formula we can transfer this information to connected bicolored graphs. 

Rather than just counting unlabeled bipartite graphs, which were counted earlier by Harary and Prins fS] 
and by Hanlon [5] , we compute their cycle index, a power series in infinitely many variables that generalizes 
both the labeled and unlabeled enumeration, and allows us to count unlabeled bipartite blocks, which has 
not been previously accomplished. 

2. The theory of species 

2.1. Introduction 

In [11], Andre Joyal introduced the notion of "species of structures", which place the classical idea of a 
"combinatorial class" (e.g., trees or permutations) in a categorical setting. A species is a functor from the 
category FinBij of finite sets with bijections to the category FinSet of finite sets with set mapqj- We write 
F[A] for the image of the set A under the species F; for example, if F is the species of bipartite graphs then 
F[A] is the set of bipartite graphs whose vertices are the elements of the set A. We refer the reader to jH 
§1.2] for details; we give here only a brief summary of the facts that we will need from the theory of species. 

The crux of species-theoretic combinatorial analysis is to focus on "labeled" combinatorial objects while 
paying careful attention to the effects of maps acting on the label set — in particular, the actions of the 
permutation groups &a acting on each set A. (We will write S„ to denote the symmetric group acting on 
the set [n] ■■= 1, 2, . . . , n, since in general only the cardinality of the label set will matter.) Treating a species 
as a functor FinBij — >■ FinSet rather than simply as a function on the class of finite sets allows us to keep 
track of this permutation information. We note that the transport F[a] of a permutation a G &a under a 
species i^ is a permutation of the set F[A] of i^-structures labeled by A. An "unlabeled" i^-structure on A 
may be understood as an orbit of i^-structures under the action of &a- 

Classical enumerative methods frequently use the algebra of generating functions, which record the num- 
ber of structures of a given size in a given class as the coefficients of a formal power series. To achieve the 
same goal in species-theoretic analysis, we define an analogous algebraic object which records information 
related to the action of the permutation groups. This object is "cycle index" of the species, a symmetric 
function defined as a sum in the power sum basis Pi = X)? ^)- (^^ some accounts of the theory the pi are 
taken simply as independent indeterminates.) 



Definition 2.1. For a species F, define its cycle index series to be the symmetric function 

z^(pi,p2,...):=E^(E MfMK>?-.-) = E^(E M^ww) 

n>0 ' ^rT6©„ ^ n>0 ' ^rre©,, ^ 

where fix(_F[(7]) :— |{s £ F[n] : F[cr]{s) = s}\, ui is the number of i-cycles of ct, and pa — p\^p'^2 



(2) 



^The use of FinSet instead of FinBij for the target category is necessary for technical reasons related to quotients. 



It is easily verified that fix(F[CT]) depends only on the cycle type of a. The cycle types of permutations 
a £ &n are in natural bijective correspondence with integer partitions X \- n (that is, weakly decreasing 
sequences (Ai, A2, . . .) such that ^^ Xi = n), and the number of permutations in (5„ of cycle type A is n\/z\, 
where if A has h parts equal to i for each i, then z\ is I'^^i! 2'^Z2! • • ■ ■ Thus we may write the sum over 
permutations in equation l^) as a sum over partitions: 

— — Z\ 

n>OX\-n ^ 

where fix_F'[A] = fix_F[cr] and p\ = Pa for any permutation a of cycle type A. 

The cycle index series Zp of a species F captures enough of its structure that we may recover both 
labeled and unlabeled enumerations. 

Theorem 2.2. The ordinary generating function F{x) for unlabeled F-structures is given by 

F{x)^Zp{x,x^,x\...). (4) 

A proof may be found in IT, §1.2]. 

Moreover, the algebra of cycle indices directly mirrors the combinatorial calculus of species. Addition, 
multiplication, and composition of species have natural combinatorial interpretations and correspond directly 
to addition, multiplication, and plethystic composition of their associated cycle indices. This last result is 
of particular interest: 

Definition 2.3. For two species F and G with G[0] ~ 0, define their composition to be the species F o G 
given by {F o G)[A] = ]\^^p(^a)^F[k] x JlseTr ^[B]) where P{A) is the set of partitions of A. 

In other words, the composition F o G is the species of F-structures of collections of G-structures. 

Definition 2.4. Let / and g be cycle indices. Then the plethysm f o g is the cycle index 

/°.9 = f{9{Pl,P2,P3,---),9{P2,P4,P6, ■■■),■■■), (5) 

where /(a, b, . . .) denotes the cycle index / with a substituted for pi, b substituted for p2, and so on. 

This is the same as the definition of plethysm of symmetric functions (see, e.g., Stanley [HI p. 447]). 
Plethysm of cycle indices then corresponds exactly to species composition: 

Theorem 2.5. For species F and G with G[0] = 0, the cycle index of their plethysm is 

ZpoG = ZpoZ(. (6) 

where o in the right-hand side is as in equation ([Sl). 

Many combinatorial structures admit natural descriptions as compositions of species. For example, every 
graph admits a unique decomposition as a (possibly empty) set of (nonempty) connected graphs, so we have 
the species identity S = £ o S'^ where £ is the species of sets, S the species of graphs, and S*^ is the species 
of connected graphs. 

The theory of species may be extended to virtual species^ which are formal differences of species. All of 
the operations for species that we have discussed extend in a straightforward way to virtual species. We 
refer the reader to [T] §2.5] for details. In particular, two virtual species F and G are compositional inverses 
ii F o G — X (or equivalently, G o F — X) where X is the species of singletons, defined by X[A] = {A} if 
\A\ = 1, and X[A] — otherwise. We write F^^^^ for the compositional inverse of F if it exists. 

2.2. F-species and quotient species 

In classical enumerative combinatorics, Burnside's lemma (also known as the Cauchy-Frobcnius lemma) is a 
powerful tool for incorporating group actions into enumerative computations. We now present an extension 
of this technique to the species-theoretic context. 



Definition 2.6. For F a finite group, a T-species F is a combinatorial species F togetiier with an action of 
F on i^-structurcs which commutes witli isomorphisms of those structures. 

For a motivating example, consider the species fcCS of fc-colored graphs; the action of 6^ on the colors 
commutes with relabelings of graphs, so fcCS is a ©^-species with respect to this action. 
From a F-species, we can construct a quotient under the action of F: 

Definition 2.7. For F a F-species, define ^/r, the quotient species of F under the action of F, to be the 
species of F-orbits of F-structures. 

A brief exposition of the notion of quotient species may be found in [T, §3.6], and a more thorough 
exposition in [5]. 

Just as with classical species, we may associate a cycle index to a F-species, following Henderson [TO] . 

Definition 2.8. For a F-species F, we define the F-cycle index Zp-. for each 7 £ F, let 

^^(t) = E^ E fix(7-i^MW (7) 

n>0 ' ae&n 

with Pa- as in equation ([4) . 

We will call such an object (formally a map from F to the ring Q[[pi,p2, • ■ •]] of symmetric functions 
with rational coefficients in the p-basis) a T-cycle index even when it is not explicitly the F-cyclc index 
of a F-species. So the coefficients in the power series count the fixed points of the combined action of a 
permutation and the group element 7. Note that, in particular, the classical ("ordinary") cycle index may 
be recovered as Zp — Zp{e) for any F-species F. 

The algebraic relationships between ordinary species and their cycle indices generally extend to the F- 
species context. The actions on cycle indices of F-species addition and multiplication are exactly as in the 
ordinary species case considered F-componentwise. The action of composition, which in ordinary species 
corresponds to plethysm of cycle indices, can also be extended: 

Definition 2.9. For two F-species F and G, define their composition to be the F-species FoG with structures 
given by {FoG)[A] = n7rGP(A)(^['''] ^ n_BG7r ^I^l) where P{A) is the set of partitions oi A and where 7 e F 
acts on a (F o G)-structure by acting on the F-structure and the G-structures independently. 



A formula similar to that of Theorem |2.5| requires a definition of the plethysm of F-symmetric functions, 
here taken from Henderson [TOl §3]. 

Definition 2.10. For two F-cycle indices / and g, their plethysm, f o g is a F-cycle index defined by 

(/ ° g)h) = f{j){9{l){Pl,P2,P3, ■ ■ .),5(7^)(P2,P4,P6, ■■■),■■ •)■ (8) 

This definition of F-cyclc index plethysm is then indeed the correct operation to pair with the composition 
of F-species: 

Theorem 2.11 (Theorem 3.1, [TU]). If A and B are T-species and B{0) = 0, then 

Z^^s^Z^oZ^. (9) 

Recall from equation ([2) that, to compute the cycle index of a species, we need to enumerate the fixed 
points of each a G ©„. To count fixed points in the quotient species ^/r we need to count the fixed F-orbits of 
a in F under commuting actions of S„ and F (that is, under an (©„ x F)-action). This may be accomplished 
by the following generalization of Burnside's lemma fl4 . (A more general result appears in ^, Theorem 

4.2b].) 

Lemma 2.12. // F and A are finite groups and S is a set with a (T x A)-action, then for any (5 G A the 
number of T -orbits fixed by 5 is -^ 'Yl-iev ^^(7; ^) ■ 

Lemma l2.12l yields a formula for the cycle index of a quotient species in terms of the F-cycle index. 



Theorem 2.13. For a T -species F, the ordinary cycle index of the quotient species ^/v is given by 

Zf/v = ^ E zlii) = ^ E ^(7 ■ ^WK- (10) 

' ' 7er ' ' ri>o 

CTee„ 

7Gr 

We will use the notation Z^ for -^ '^■y^r ^'pil)- 

Note that this same result on cycle indices is implicit in [21 §2.2.3]. With it, we can compute explicit 
enumerative data for a quotient species using cycle-index information of the original F-species with respect 
to the group action. 

We note that when F is a symmetric group ©„, as in this case, we may represent the F-cycle index as a 
symmetric function in two sets of variables which is homogeneous of degree n in the second set of variables; 
with this approach, F-cyclc index plethysm corresponds to the operation of "inner plethysm in y" studied 
by Travis [17]. 

3. The species of bipartite blocks 

3.1. Introduction 

Definition 3.1. A hicolored graph is a graph G each vertex of which has been assigned one of two colors 
(here, black and white) such that each edge connects vertices of different colors. A bipartite graph (sometimes 
called bicolorable) is a graph G which admits such a coloring. 

There is an extensive literature about bicolored and bipartite graphs, including enumerative results for 
bicolored graphs [7], bipartite graphs both allowing 5j and prohibiting [8_ isolated points, and bipartite 
blocks [S]- However, the enumeration of bipartite blocks has been accomplished previously only in the 
labeled case. By considering the problem in light of the theory of F-species, we develop a more systematic 
understanding of the structural relationships between these various classes of graphs, which allows us, in 
particular, to enumerate all of them in their unlabeled forms. 

Throughout this chapter, we denote by 23C the species of bicolored graphs and by "BCP the species of 
bipartite graphs. The prefix C will indicate the connected analogue of such a species, so CCBIP is the species 
of connected bipartite graphs. 

We are motivated by the graph-theoretic fact that each connected bipartite graph has exactly two bicol- 
orings, and may be identified with an orbit of connected bicolored graphs under the action of ©2 where the 
nontrivial element r reverses all vertex colors. We will hereafter treat all the various speci es o f bicolored 



graphs as S2-species with respect to this action and use the theory developed in Section |2.2| to pass to 
bipartite graphs. 

3.2. Bicolored graphs 

We begin our investigation by directly computing the (32-cycle index for the species 'BC of bicolored graphs 
with the color-reversing ©2-action described previously. We will then use various methods from the species 
algebra of Section |2 to pass to other related species. To compute the 62-cycle index Z^g we compute 
separately Z®g(e) and Z®g(r). 

3.2.1. Computing Z||(e) 

We construct the cycle index for the species 23C of bicolored graphs in the classical way, which in light of 
our (32-action will give Z^^{e). 

For each n > and each permutation tt G 6„, we must count bicolored graphs on [n] for which tt is 
a color-preserving automorphism. To simplify some future calculations, we omit empty graphs and define 
23C[0] = 0. We note that the number of such graphs in fact depends only on the cycle type A h n of the 
permutation tt, so we can use the cycle index formula in equation ([3|) interpreted as a F-cycle index identity. 



Fix some n > and let X h n. We wish to count bicolored graphs for which a chosen permutation tt 
of cycle type A is a color-preserving automorphism. Each cycle of the permutation must correspond to a 
monochromatic subset of the vertices, so we may construct graphs by drawing bicolored edges into a given 
colored vertex set. If we draw some particular bicolored edge, we must also draw every other edge in its 
orbit under tt if tt is to be an automorphism of the graph. Moreover, every bicolored graph for which tt is 
an automorphism may be constructed in this way Therefore, we direct our attention first to counting these 
edge orbits for a fixed coloring; we will then count colorings with respect to these results to get our total 
cycle index. 

Consider an edge connecting two cycles of lengths m and n; the length of its orbit under the permutation 
is lcm(?Ti, n), so the number of such orbits of edges between these two cycles is mn/ lcm(TO, n) = gcd(m, n). 
For an example in the case m = 4, ra = 2, see Fig. y. The number of orbits for a fixed coloring is then 
^ gcd(TO, n) where the sum is over the multisets of all cycle lengths m of white cycles and n of black cycles 
in the permutation tt. We may then construct any possible graph fixed by our permutation by making a 
choice of a subset of these cycles to fill with edges, so the total number of such graphs is Y[ 2S'='^(™^") for a 
fixed coloring. 



.---Q 




Figure 2: An edge e (solid) between two cycles of lengths 4 and 2 in a permutation and that edge's orbit 
(dashed) 

We now turn our attention to the possible colorings of the graph which are compatible with a permutation 
of specified cycle type A. We split our partition into two subpartitions, writing X = fiU v, where partitions 
are treated as multisets and U is the multiset union, and fi represents the white cycles and v the black. Then 
the total number of graphs fixed by such a permutation with a specified decomposition is 



fix(M,i.)=n2^''^''' 



j) 



(11) 






where the product is over the elements of /i and A taken as multisets. Suppose that the multiplicities of the 
part i in the partitions A, /x, and u are h, rrii, and n^, respectively. Then the U i-cycles of a permutation of 
cycle type A can be colored so that rrii are white and rii are black in Zi!/(mi!ni!) ways. So in all there are 
Yii h^-fi'mi^- ni\) — z\/(Zfj_z,y) colorings associated with fi and v, and 



fix(A) 



zx 



fix(/i, u) 



jei' 



fiUu—X ^ i^f2 



Thus we have a formula for Z,^^{e): 
Theorem 3.2. 



zsm=Y: E ^n^^"'"^ 



M,l^j) 



„>0 fJ.,!^ t^ ^ jj 



(12) 



Explicit formulas for the generating function for unlabeled bicolored graphs were obtained by Harary 
[7] using conventional Polya-theoretic methods. Conceptually, this enumeration in fact largely mirrors our 
own. Harary uses the classical cycle index of the line groupu of the complete bicolored graph of which any 
given bicolored graph is a spanning subgraph. He then enumerates orbits of edges under these groups using 
the Polya enumeration theorem. This is clearly analogous to our procedure, which enumerates the orbits of 
edges under each specific permutation of vertices. 

3.2.2. Calculating Z^I{t) 

Recall that the nontrivial element of r e 62 acts on bicolored graphs by reversing all colors. 

We again consider the cycles in the vertex set [n] induced by a permutation tt S ©« and use the partition 
A corresponding to the cycle type of n for bookkeeping. We then wish to count bicolored graphs on [n] 
for which t • tt is an automorphism, which is to say that tt itself is a color-reversing automorphism. The 
number of bicolored graphs for which tt is a color-reversing automorphism depends only on the cycle type A. 
Each cycle of vertices must be color-alternating and hence of even length, so the partition A must have only 
even parts. Once this condition is satisfied, edges may be drawn either within a single cycle or between two 
cycles, and as before if we draw in any edge we must draw in its entire orbit under tt (since tt is to be an 
automorphism of the underlying graph) . Moreover, all graphs for which tt is a color-reversing automorphism 
with a fixed coloring may be constructed in this way, so it suffices to count such edge orbits and then consider 
how colorings may be assigned. 

We first determine the number of orbits of edges within a cycle of length 2n; we hereafter describe such 
a cycle as having semilength n. There are exactly n^ possible white-black edges in such a cycle. If n is even, 
then every edge lies in an orbit of size 2n, so there are n^/{2n) — n/2 orbits of edges. If n is odd, there are 
n edges joining diametrically opposed vertices, which have oppositive colors. These n edges are all in the 
same orbit. (See Fig. |3a| for an illustration of these edges.) The remaining n^ — n edges are in orbits of size 



2n, so there are (n — n)/{2n) — [n — l)/2 of these orbits. (See Fig. |3b| for an illustration of these edges.) 



Thus the total number of orbits for n odd is (n + l)/2. In either case, the number of orbits is \n/2~\ . 

..q"' ■■•p. .a- 

\ / 

^^^D # 



/ \ 
/ \ 



"■tj.... ,..b-" "■t>.:.:.:::-£)-" 

(a) A diameter d {I = 3) (b) A non-diameter e {I = 1) 

Figure 3: Both types of intra-cycle edges and their orbits on a typical color-alternating 6-cycle 

Now consider an edge drawn between two cycles of semilengths m and n. The total number of possible 
white-black edges is 2mn, each of which has an orbit length of lcm(2?Ti, 2n) = 21cm(m, n). Hence, the total 
number of orbits is 2mn/{2lcTa{m,n)) = gcd(m, n). 

All together, then, the number of orbits for a fixed coloring of a permutation of cycle type 2 A (denoting 
the partition obtained by doubhng every part of A) is J2i \^i/'^^ + X]i<j gcd(Ai, Aj). All valid bicolored 
graphs for a fixed coloring for which tt is a color-preserving automorphism may be obtained uniquely by 
making some choice of a subset of this collection of orbits, just as in Section l3. 2. ll . Thus, the total number 



^The line group of a graph is the group of permutations of edges induced by permutations of vertices. 
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Figure 4: An edge e and its orbit between color-alternating cycles of semilengths 2 and 1 



of possible graphs for a given vertex coloring is 



i i<j 



(13) 



which we note is independent of the choice of coloring. For a partition 2 A with l{\) cycles, there are then 
2H-*') colorings compatible with our requirement that each cycle is color-alternating, which we multiply by 
(I13p to obtain the total number of graphs for all permutations tt with cycle type 2A. Therefore we have a 
formula for Z®|(t): 



Theorem 3.3. 



^se 



M-E E 



n>0 Ahr 
n even 



./2 



2i(A)P2A 
Z2\ 



J|2rA,/2l-Q2' 



gcd(Ai,Aj) 



(14) 



i<j 



3.3. Connected bicolored graphs 

As noted in the introduction to this section, we may pass from bicolored to bipartite graphs by taking a 
quotient under the color-reversing action of ©2 only in the connected case. Thus, we must first pass from the 
species 236 to the species CSC of connected bicolored graphs. It is a standard principle of graph enumeration 
that a graph may be decomposed uniquely into (and thus species-theoretically identified with) the set of its 
connected components. This same relationship holds in the case of bicolored graphs. Thus, the species CBC 
of nonempty bicolored graphs is the composition of the species CSC of nonempty connected bicolored graphs 
into the species £+ = £ — 1 of nonempty sets: 



SC = £+oCSC. 



(15) 



Reversing the colors of a bicolored graph is done simply by reversing the colors of each of its connected 
components independently; thus, once we trivially extend the species £"'' to an 62-species by applying the 
trivial action, equation (J15l ) holds as an identity of (32-species for the color-reversing (32-action described 
previously. 

To use the decomposition in equation (jlSl ) to derive the 62-cycle index for CBC, we must invert the 



62-species composition into £^. We write O := (£ 



+ ^(-i> 



to denote the virtual species that is the inverse 



of £+ with respect to composition of species, following the notation of y^. We can derive from pj §2.5, 
equation (58c)] that its cycle index is 



Zr 



E^iogd 



k>l 



k 



-Pk) 



(16) 



where /z is the Mobius function. We can then rewrite equation (J15r ) as 

C!BC = r2o'BC. 



From Theorem |2. Ill we get a formula for the (32-cycle index for CSC: 



Theorem 3.4. 



^ese — ^n ° ^se- i^'^) 



Note that we could have avoided the use of virtual species by performing the inversion at the level of 
cycle indices. 

3.4. Bipartite graphs 

As we previously observed, connected bipartite graphs are naturally identified with orbits of connected 
bicolored graphs under the color-reversing action of 62- Thus, 



By application of Theorem |2.13| . we can then directly compute the cycle index of CBT in terms of previous 
results. 

Theorem 3.5. 



Since a bipartite graph is a set of connected bipartite graphs, we have HT = £ o 623 T, and this gives the 
formula for the cycle index for bipartite graphs. 

Theorem 3.6. 

Z^y — Z^O Zg3y. (19) 



Theorem |3.6| allows us to compute the number of unlabeled bipartite graphs with n vertices. However, 
to do this calculation, we don't need the entire cycle index; we can do all of the calculations using ordinary 
generating functions. Specifically, let 

fe{x)^l + Z'il{e){x,x\x\...) 

fr{x)^l + Z^l{T){x,x\x\...) 

ge{x) = Zg|g(e)(a;,a;^a;^,...) 
9t{x) — Zq^q[t){x,x ,x ,...) 

c{x) = Ze^vix) = Z(o,^,j,{x, x'^,x^, . . . ) 

b{x) = Z'B-j>{x) = Z,_g.jy{x, x'^,x^, . . .). 

Then c{x) is the ordinary generating function for connected bipartite graphs and b{x) is the ordinary gener- 
ating function for bipartite graphs. We have formulas for feix) and frix) as sums over partitions. 






gcd(/Ji,i/j) 



ri=0 p. I' ^ " i,j 



n even Ahn/2 " i i<j 



and ge{x) and grix) are related to fe{x) and frix) by 






«-)=»pi:^i^+e 



_.. , ^ 2k 

k=0 k=l 



which may be inverted to give 

5e(-)=f:^l0g/e(x'=) 

9Ax) = f2t^^§±^logfA^''^') + f:^logUx'% (20) 

fe=0 fe=l 

where /x is the Mobius function. FinaUy, c{x) — ^{ge{x) + grix)) and b{x) — exp(^^^ c(a;'^)/fc). 

These calculations are essentially the same as Hanlon's [5), though he does not have our equation (|2C 
Unlabeled bipartite graphs were first counted, using a different approach, by Harary and Prins [5]. 

In order to count bipartite blocks, which we accomplish in the next section, we do need the entire cycle 
index. 

3.5. Nonseparable graphs 

We now turn our attention to the notions of block decomposition and nonseparable graphs. A graph is said 
to be nonseparable if it is vertex- 2-connected (that is, if there exists no vertex whose removal disconnects the 
graph); every connected graph then has a canonical "decomposition'lfl into maximal nonseparable subgraphs, 
often shortened to blocks. In the spirit of our previous notation, we we will denote by DsfST the species of 
nonseparable bipartite graphs, our object of study. 

The basic principles of block enumeration in terms of automorphisms and cycle indices of permutation 
groups were first identified and exploited by Robinson \i^. In [1, §4.2], a theory relating a species B of 
nonseparable graphs to the species Cb of connected graphs whose blocks are in B is developed using similar 
principles. 

We extract two particular results, appearing as |1] equations 4.2.27 and 4.2.26a]. We note that the 
derivative F' of a species F 1, pp. 47-49] is defined by F'[A\ — F[A U {*}], where * is not in A, and its 
cycle index satisfies Zp' = dZp/dpi. The pointing F* of F [T, §2.1] is XF'. 

Lemma 3.7. Let B be a species of nonseparable graphs and let C denote the species of connected graphs 
whose blocks are in B. Then 

B^cfc'^-^A +XB' ^X (21a) 

and 

It is apparent that the class of nonseparable bipartite graphs is itself exactly the class of blocks that occur 
in block decompositions of connected bipartite graphs. Accordingly, we can apply Lemma |3/7| immediately: 



Theorem 3.8. 



by equation (|21al ) and 



'N'By = G'ByOs'By'^-^^) + x ■ j^'BT' - x, (22a) 



^^^' = "(^^^^1 (22b) 



by equation (|21 



We have already calculated the cycle index for the species C23CP, so the calculation of the cycle index of 
J^'By is now simply a matter of algebraic expansion. 

A generating function for labeled bipartite blocks was given by Harary and Robinson [9], where their 
analogue of equation |22l ) for the labeled exponential generating function for blocks comes from ^4]. However, 
we could locate no corresponding unlabeled enumeration in the literature. The numbers of labeled and 
unlabeled nonseparable bipartite graphs for 71 < 10 as calculated using our method are given in Table [ij. 



^Note that this decomposition does not actually partition the vertices, since many blocks may share a single cut-point, a 
detail which significantly complicates but does not entirely preclude species-theoretic analysis. 
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A. Cycle indices of compositional inverse species 

In Section Is. a . our results included two references to the compositional inverse CBT*'^^' of the species CSCP*. 
Although we have not explored computational methods in depth here, the question of how to compute the 
cycle index of the compositional inverse of a specified species efficiently is worth some consideration. Several 
methods are available, including one developed in 1, 4.2.19] as part of the proof that arbitrary species have 
compositional inverses, but our preferred method is one of iterated substitution. 

Suppose that ^ is a species (with known cycle index) of the form X + 'i>2 + "^3 + ■ ■ ■ where ^i is the 
restriction of ^ to structures on sets of cardinality i and that $ is the compositional inverse of ^P. Then 
* o $ = X, so 

X = ^o$ = $ + *2(*) + 'i'si^) + ■■■ 

also. Thus 

$ = X-*2($) -*3($) -••■■ (23) 

This recursive equation is the key to our computational method. To compute the cycle index of $ to degree 
2, we begin with the approximation $ « X and then substitute it into the first two terms of equation (|23l ): 
$ « X — '^2{X) and thus Z^ « Z^^ — Z^ o Z-^. All terms of degree up to two in this approximation will be 
correct. To compute the cycle index of $ to degree 3, we then take this new approximation $ w X — 'i>2{X) 
and substitute it into the first three terms of equation (J23r ) . This process can be iterated as many times as are 
needed; to determine all terms of degree up to n correctly, we need only iterate n times. With appropriate 
optimizations (in particular, truncations), this method can run very quickly on a personal computer to 
reasonably high degrees; we were able to compute Z (_i) to degree sixteen in thirteen seconds. 

B. Numerical results 

With the tools developed in Section |3|, we can calculate the cycle indices of the species AfBV of nonseparable 
bipartite graphs to any finite degree we choose using computational methods. This result can then be used to 
enumerate unlabeled bipartite blocks. We have done so here using Sage [16] and code listed in Appendix O 
The resulting values appear in Table [ij. 
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Table 1: Enumerative data for unlabeled bipartite blocks with n hedra 



n 


Unlabeled 


1 


1 


2 


1 


3 





4 


1 


5 


1 


6 


5 


7 


8 


8 


42 


9 


146 


10 


956 


11 


6643 


12 


65921 


13 


818448 


14 


13442572 


15 


287665498 


16 


8099980771 


17 


300760170216 


18 


14791653463768 


19 


967055338887805 


20 


84368806391412395 


21 


9855854129239183783 


22 


1546801291978378704267 


23 


327092325302250220001201 


24 


93454432085788531687319514 
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C. Code listing 

The functional equation (|22f ) characterizes the cycle index of the species TsfBCP of bipartite blocks. In this 
section we have used the the computer algebra system Sage [T^ to adapt the theory into practical algorithms 
for computing the actual numbers of such structures. Python/Sage code to compute the coefficients of the 
ordinary generating function 'N'B7'{x) of unlabeled bipartite blocks explicitly follows in listing [l|. 

Listing 1: Sage code to compute numbers of bipartite blocks (bpblocks .sage) 

1 # Import needed code 

from sage . combinat . species . stream import Stream, _integer s_f rom #Infinite generator for lazy 
power series 
3 from sage . combinat . species . generating_series import CyclelndexSer iesRing 

5 # Set up helper variables 

# We'll work with these cycle indices a lot 

7 X = species . SingletonSpecies (). cycle_index_ser ies () 

E = species . SetSpecies (). cycl6_index_ser ies () 
9 

CIS = CyclelndexSeriesRing (QQ) #The ring of cycle index series with rational coefficients 

11 p = SFAPower(QQ) #The ring of symmetric functions (power-sum basis) with rational coefficients 

13 Z2 = Symmetr icGroup (2) #The group of order 2 

e = Z2 . identity () 
15 t = Z2 . gen () 

17 # Set up some helper functions 

# Compute the derivative of a specified cycle index 
19 def ci_deriv ( f ) : 

# Helper function to build a CIS corresponding to the ith term of the derivative of / 
21 def termbuilder ( i ) : 

# The first i coefficients are 0, the (i + l)th is the derivative of the 
23 # corresponding coefficient of /, and those afterwards are (filled 

# automatically by the CIS constructor) . 

25 helperlist = [0]*i + [f . coefficient (i+1) . derivat ive_with_resp6ct_to_pl () , 0] 

# Now we let the CyclelndexSeriesRing constructor build a suitable CyclelndexSeries 
27 return CIS (helperlist ) 

29 # The sum_generator constructor then builds a CyclelndexSeries which sums 

# all the terms we constructed above . 

31 return CIS . sum_generator ( termbuilder ( i) for i in _integers_f rom (0) ) 

33 # Compute the pointing of a specified cycle index 

ci_pointed = lambda f: ci_deriv(f) * X 
35 

# Compute the compositional inverse of a specified cycle index 
37 # (Algorithm from Appendix 0) 

def ci_compinv ( f ) : 
39 # We must have that 3^0 = 

assert f . coefficient (0) == X . coefficient (0) 
41 # and that 'Ji = X 

assert f . coefficient (1) == X . coefficient (1) 
43 # in order to apply this algorithm. 

45 # Following the result in Appendix Q, if the conditions above are satisfied, 

# we have ?<-!> = X- (J - X) o (J^-l)). 

47 # The Sage cycle index define method handles this recursive definition. 

result = CISO 
49 result . def ine (X - (f - X) . compose (result ) ) 

51 return result 

53 # Compute the multiplicative inverse of a specified cycle index 

# (Algorithm from JT], §2.5 exercise 7]) 
55 def ci_multinv ( f ) : 

# The result of Q requires that Jo 5^ . 
57 assert f . coefficient (0) != 
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59 # If this condition holds, we have that 3""^ = ^°^„{— ^o)* ■ {^+)*. just as for 

# ordinary formal power series. 
61 return CIS.sum_generator(f.coefficient(0)"(-i-l) * (f. coefficient (0) - f)''i for i in 

_integers_f rom (0) ) 

63 # Compute a cycle index divided by X, if possible 

# (Helper for a later calculation where J~^ is not defined but t^ is) 
65 def ci_xdiv ( f ) : 

def pl_dropper ( part ) : 
67 assert 1 in part 

return p(part[:-l]) 
69 

def pl_dropper_sf ( sf ): 
71 assert sf in SFAPower(QQ) 

return SFAPower (QQ ) . _apply_module_morphism (sf , pl_dropper) 
73 

termbuilder = lambda i: ClS([0]*i + [pl_dropper_sf (f . coefficient (i+1) ) , 0]) 
75 return CIS . sum_generator ( termbuilder ( i) for i in _integers_f rom (0) ) 

77 # Compute the cycle index for the combinatorial logarithm Q 

# (Algorithm from II12L Prop. 6] ; somewhat faster applying the more general inversion algorithm above to 

compute E^^^^f) 
79 

# First we define a helper function to compute the terms of Z^-. . 
81 def omegaterm (n) : 

assert n>=0 
83 assert n in ZZ 

85 if n == 0: 

return 
87 elif n == 1: 

return p [1] 
89 else: 

return 1/n * ((-l)*(n-l) * p[l]**n - sum(d * 

p ( [Integer (n/d) ]). plethysm ( omegaterm (d) ) for d in di visor s (n) [: -1] ) ) 
91 

# Then we convert this to a Python generator for the benefit of the Sage 
93 # CyclelndexSeries constructor. 

def omegagen () : 
95 for n in _integers_f rom (0) : 

yield omegaterm (n) 
97 

# Then we feed this generator into our CyclelndexSeriesRing. 
99 Omega = CIS ( omegagen () ) 

101 # Helper functions for working with partitions 

# The union of two partitions is the partition corresponding to their multiset union 
103 def union ( mu , nu) : 

return Partition ( sorted (mu . to_list () + nu . to_list () , reverse=true ) ) 
105 

# For a partition ft = [fJ-i,fJ.2, . . .] h m and a natural n, we define n ■ /.t := [nij,i,nfj.2, . . .] h nm. 

107 # That is, /; ■ // is the partition obtained from /i by multiplying the multiplicity of each part of /.( by n. 

def partmult ( mu , n ) : 
109 return Partition ( [part * n for part in mu . to_list () ] ) 

111 # Define the cycle index for 256 

# First we compute the number of graphs fixed un der e by permutations of cycle 
113 # types fj, and u in accordance with equation dill ) . 

def ef ixedbcgraphs ( mu , nu ): 
115 return 2** ( sum ( [gcd (i , j) for i in mu for j in nu] ) ) 

117 # Then we build a generator for the terms of the cycle index Z^g(e). 

def egen () : 
119 yield p(0) 

for n in _integers_f rom (1) : 
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121 yieldsumCp (union (palr[0], pair[l]))/(pair[0]. ant () * pair[l]. aut () ) * 

ef ixedbcgraphs (pair [0] , pair[l]) for pair in Part it ionTuples (n , 2)) 

123 # Then we compute the number of grap hs fixed under r by a permutation of cycle type 

# /J in accordance with equation dlSl ) . 
125 def tf ixedbcgraphs ( mu ): 

return 2**(len(mu) + sum ( [ integer_ceil (p/2) for p in mu] ) + sum ( [gcd (mu [i] , mu[j]) for 
i in range(0, len(mu)) for j in range(i+l, len(mu))])) 
127 

# Then we build a generator for the terms of the cycle index Z^^[t). 
129 def tgenO : 

yield p(0) 
131 for n in _integers_f rom (1) : 

yield p(0) 
133 yield sum (tf ixedbcgraphs ( mu ) * p ( part mult (mu , 2))/partmult (mu , 2).aut() for mu 

in Part it ions (n) ) 

135 # Finally, we use the CyclelndexSeriesRing constructor (twice) to define Z,'^'?, as a dictionary. 

BC = {e: CIS(egen()), t: CIS(tgen())} 
137 

# Define the cycle index for CBS 

139 # The S2-cycle index plethysm Z^^ o Z^?, requires 

# the more complex calculation in Definition |2. lOl . which we compute termwise. 
141 def CBCtermmap (term) : 

if term == 0: 
143 return CIS (0) 

145 # We use the r'' term stretched by p for each part p in our partition 

termbuilder = lambda part: prod (BC [t**p] . stretch (p) for p in part) 

147 

# We then return a CyclelndexSeries summing the new terms with their original coefficients 

149 return sum ( coeff *termbuilder (part ) for part , coef f in term) 

151 # The e term of -^g^g is a classical plethysm. 

# The T term uses our more sophisticated calculation above. 

153 CBC = {e: Omega . compose (BC [e] ) , t: CIS . sum_generator ( CBCtermmap ( Omega . coef ficient ( i) ) for 
i in _integers_f rom (0) ) } 

155 # Define the cycle indices for C'&y and SCP as in Theorem |3 

CBP = (CBC[e] + CBC[t]) * (1/2) 
157 

# Define the cycle index for 21" as in Theorem |3.6| . 
159 BP = E . compose (CBP) 

161 # Compute the cycle index for ?^2^T as in Theorem |3.8| . 

CBP_pointed_compinv = ci_compinv ( ci_pointed (CBP) ) 
163 NBP = CBP . compose (CBP_pointed_compinv) + X * 

Omega. compose (ci_multinv (ci_xdiv (CBP_pointed_compinv) ) -1) 

165 # Compute the first k coefficients of Ji1i'J'{x) 

# (Note that this list is zero-indexed!) 

167 NBP . isotype_generating_series () . coefficients (k) 
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